Stimulated Raman adiabatic passage in an open quantum system: Master equation 

approach 



M. Scala, 1 B. Militcllo, 1 A. Messina, 1 and N. V. Vitanov 2 ' 3 

1 Dipartimento di Scienze Fisiche ed Astronomiche dell'Universita di Palermo, Via Archirafi 36, 90123 Palermo, Italy 
2 Department of Physics, Sofia University, James Bourchier 5 blvd, 1164 Sofia, Bulgaria 
'^Institute of Solid State Physics, Bulgarian Academy of Sciences, Tsarigradsko chaussee 72, 1784 Sofia, Bulgaria 

A master equation approach to the study of environmental effects in the adiabatic population 
transfer in three-state systems is presented. A systematic comparison with the non-Hermitian 
Hamiltonian approach [N. V. Vitanov and S. Stenholm, Phys. Rev. A 56, 1463 (1997)] shows that 
in the weak coupling limit the two treatments lead to essentially the same results. Instead, in the 
strong damping limit the predictions are quite different: in particular the counterintuitive sequences 
in the STIRAP scheme turn out to be much more efficient than expected before. This point is 
, explained in terms of quantum Zeno dynamics. 
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I. INTRODUCTION 



Stimulated Raman adiabatic passage (STIRAP) is a powerful technique for coherent population transfer in a three- 
state chainwise connected system 1-2-3 [H-@. In the most common A linkage pattern, states |1) and |3) are ground 
*^ or mctastablc levels, while the intermediate state |2) is a decaying excited electronic level. The unique advantage of 
II STIRAP over other population transfer techniques is that in the adiabatic limit this intermediate state |2) does not 
receive even transient population during the transition |1) — > |3). This feature derives from the fact that STIRAP 
proceeds via a dark state, which is a superposition of states |1) and |3) only. A two-photon resonance between states 
^ ' |1) and 1 3) ensures the emergence of such an eigenstate of the Hamiltonian. A counterintuitive pulse sequence, Stokes- 
O" 1 pump (with the Stokes driving the 2-3 transition and the pump driving the 1-2 transition), aligns initially state |1) 
with the dark state, Finally, adiabatic evolution, which is enforced by selecting sufficiently large pulse areas, ensures 
that the three-level system remains in the dark state at all times until it aligns with the target state |3) in the end. 
In the adiabatic limit, the properties of the intermediate state |2), including its detuning and loss rates (e.g., 
^sjQ ' spontaneous emission within and outside the system), arc irrelevant because it is decoupled from the dynamics. 
i~H However, these factors cannot be ignored completely because first, in a real experiment the evolution is never perfectly 
^^O ' adiabatic, and second, these factors affect the adiabatic condition itself. The robustness of STIRAP against the 
intermediate- level detuning has been quantified in Q- The effects of dephasing Q and spontaneous emission within 
C*~) ' the system have also been scrutinized. 

The effect of irreversible population loss from the intermediate state |2) has been studied in [l(| wherein the decay 
■ has been introduced phenomenologically, adding an imaginary diagonal term in the lossless Hamiltonian. Nevertheless, 
because of the rapidly increasing popularity of STIRAP as a quantum control tool in dissipative environments, 
it is interesting, instructive and important to treat the problem with greater mathematical rigor, starting from a 
microscopic model, which explicitly takes into account the coupling between the system and an external environment. 
For time-independent models, there are well established techniques which allow for a description of the open 
' dynamics of the system of interest by means of master equations, which can be systematically derived from the 
system Hamiltonian [Tll.[l2l|. However, for time-dependent Hamiltonian models the derivation of the master equation 
requires more attention. A fully satisfactory and very simple theory of master equations for such systems has been 
developed by Davies in the late 70's [l3|. The main feature of this approach is that, under the hypothesis of very short 
reservoir correlation times, one obtains a time-dependent master equation describing jumps between instantaneous 
cigenstates of the system Hamiltonian. This approach has been used recently in the study of quantum logic gates 
based on the accumulation of geometric phases in adiabatic evolutions [3] ■ Similar time-dependent master equations 
have been used by other authors in problems involving quantum adiabatic evolution [TEI . [la ] . 

We emphasize that in general the microscopic derivation of a master equation for a given physical system may give 
rise to predictions which differ from the ones obtained from phenomenological models, as recently seen for example 
in the context of lossy cavity QED [l7l - f2"(jj and two-qubit dynamics [2l[ . 

In this paper, we present a microscopic model from which we derive a master equation describing the dissipative 
dynamics of a system subjected to a STIRAP scheme. In order to perform a systematic comparison between our 
model and the phenomenological model in [lj| , we move to a description of the dynamics by means of an effective non- 
Hcrmitian Hamiltonian, in this case equivalent to the master equation approach. Comparing the predictions coming 
from the two non-Hermitian Hamiltonian models, we find that, according to our effective model, the population 
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FIG. 1: (Color online) Level scheme and coupling scheme. Level 4 is very far from the other three. 



transfer is more efficient than previously expected. The discrepancy is more evident in the limit of strong damping, 
where the new effects found can be easily understood in terms of quantum Zeno dynamics. 

The paper is structured as follows. In the next section we recall the main properties of the system under scrutiny 
which has been described in [Io| and derive the master equation starting from a microscopic model of system-reservoir 
coupling. In the third section we derive the relevant effective Hamiltonian, while in the fourth section we compare 
the predictions from the effective and the phenomenological Hamiltonian models. Finally, in the last section some 
conclusive remarks are given. 



II. MASTER EQUATION 

We model the irreversible population decay from the intermediate state |2) by introducing an additional state |4), 
to which state |2) decays (see figure In the absence of such a decay the Hamiltonian describing the enlarged 
four-state system reads (with h = 1): 



H (t) = 
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with W4 < well separated from the other levels, and ui mn = uj m — ui n . The population decay of state |2) to state 
|4) is modeled by a system-bath coupling Hamiltonian involving these two states, which in the non-rotating frame is 
given by: 



#dcc = (|2) (4| + |4) (2|) ® Y,9k (ok + 4 



(2) 



In the rotating frame, the transformation to which is given by T(t) = e iwi * |1) (1| + e *("2-A)t | 2 ) (2\ + e l ^t ^ ^ ^ 
Hamiltonian of the system becomes 



and -ffdec changes accordingly: 
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H dcc = [e l ^ 2 ' A ) t |2)(4|+e- l( " 2 - A)t |4)(2| 



(4) 
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The instantaneous eigenstates of the Hamiltonian H s (t) are [l(| : 



|+) = sin sin |1) + cosy? |2) + sin ^ cos 13) , (5a) 
|0) =cos(9|l) -sin<9|3) , (5b) 
|— ) = cos ip sin 8 |1) — sin ip |2) + cos ip cos 6* |3) , (5c) 



where 



tan*(*) = £|, (6a) 
tan2^) = ^|), (6b) 

fio(t) = y / n p (t) 2 + O s (t) 2 , (6c) 

and state |4) which is left unchanged by the transformation. The corresponding eigenvalues are cj + = Slocoty), 0, 
cj_ = — Ootany? and u^. 

In the ideal case (in the absence of i?dcc)i perfect population transfer from |1) to |3) takes place. In particular, two 
different pulse sequences are possible, the intuitive and the counterintuitive sequences. In the intuitive sequence, the 
pump pulse fl p (t) precedes the Stokes pulse Fl s (t); then for nonzero single-photon detuning A it can be shown that 
in the adiabatic limit the population remains at all times in state |— ), which at t = — oo is equal to state |1) and 
at t = oo it is equal to state |3). Therefore, in the adiabatic limit there is a perfect population transfer from |1) to 
|3) [lU; this process has recently been termed b-STIRAP (because it proceeds via the bright state |— )) (23|. For the 
counterintuitive pulse sequence, the pump pulse follows the Stokes pulse. Then the population is transferred through 
the dark state |0), which again is equal to |1) at t = — oo and to |3) at t = oo. 

In the non-ideal case of decaying state |2), the predictions change. In ref. [l(| it has been shown that the instability 
of the intermediate state reduces the efficiency of the scheme. In particular it turns out that the intuitive sequence is 
more fragile than the counterintuitive one. Since these results have been obtained by means of a phcnomcnological 
Hamiltonian approach, it is natural to ask whether these results may change when a microscopic master equation 
approach is used. In the following we give the master equation describing the dissipative dynamics of the system. 

According to the general theory by Davies [l3j|, recently applied to some time-dependent systems under the 
assumption of bath correlation times much smaller than the time of variation of the Hamiltonian, the master equation 
describes jumps between the instantaneous eigenstates of the time-dependent Hamiltonian. Since the Hamiltonian 
(0| involves only states |2) and |4), and since state |2) is not involved in the dark state |0), the only jumps allowed 
are |+) |4) and |— ) ^ |4). Therefore the master equation of our system is given by 



dp 
dt 



= -i[H s (t),p] 

+ 7+ (t)(|4)(+M+)(4|-i{l+}<+l,P}) 

+ 7 _(*)(|4)HpI->(4|-HI->H,p}) 

+ 7 rW(l+)(4|p|4)(+|-i{|4)(4|,p}) 

+ 7 CTC W(h)(4|p|4)<-|-i{|4)(4|,p}) • (7) 



Here 



7+ (t) = cos 2 p(t)\g(u + )\ 2 D(u + )(N(u+) + 1) , (8a) 
7 _(i) = sin 2 <p(t)\g(£j_)\ 2 D(£j_)(N(£j_) + 1) , (8b) 

are the decay rates of the states |+) and |— ) toward state |4) , corresponding to the Bohr frequencies u>± = uj± +uj2 — ^i- 
The quantities D(uj±) and N(uj±) are respectively the densities of modes and the average numbers of photons in 
the reservoir at the relevant frequencies. The parameter g{oJ±) gives the continuum-limit system-reservoir coupling 
strengths. We note that the factors cos 2 <p and sin 2 tp in Eq. © come from the calculation of the square moduli of 
the matrix elements of the system operator |2) (4| + |4) (2| appearing in -ffdoc between states |4) and |+), and between 
|4) and |— ), respectively. Finally, the excitation rates are given by 7± BC (t) = r y±(t)N(ui±)/(N(uj±) + 1) and vanish at 
zero temperature. This will be the case in the following sections. 
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III. RATE EQUATIONS AND EFFECTIVE HAMILTONIAN 

Introducing the density matrix decomposition in terms of the instantaneous eigenstates of the Hamiltonian H s (t), 

p(t) = Y,pMm(m\, (9) 



and substituting it into the master equation ([7]), one obtains a set of rate equations, 

Pij(t) = ^A ijk i{t)p M {t) ■ 



(10) 



On the other hand, introducing the matrix of the density operator matrix elements in the time-dependent basis, i.e. 

P =[/*,(*)], (11) 

with i,j = ±,0, one can prove (see appcndix|X| that the same rate equations can be obtained with a pseudo-Liouvillc 
equation, restricted to the subspace |+), |0), |— ): 



with 



Hcfi = 



p = -i( H oS p - pHl & , 



f2o cot p — iT cos 2 p iO sin ip ip 

—iOsmp — iOcosp 

—ip iO cos ip — Qq ta,np — irsin 2 p> 



(12) 



(13) 



where, reminding that N(uj±) = at zero temperature and assuming flat reservoir spectrum jl2| . we have set 

r = |. 9 (^ + )| 2 J D(^ + ) = |^_)| 2 ^). 

The comparison of this Hamiltonian with the phenomcnological Hamiltonian of Eq. (4) of [1^, i.e.: 



where 



phcn 



f2o cot tp — iT cos 2 p i6 sin 93 iip + ^rsin2i^ 
—id sin ip — iO cos p> 

— iip + iT sin2ip iO cos p —VIq tm\ p — iT siv? p 



(14) 



(15) 



shows that the only difference is in terms (1,3) and (3,1). The two Hamiltonians (fT3")) and ([T3|) are written down in 
the basis |+), |0), I-) (not |1), |2), |3)). 

More details of the derivation of the effective Hamiltonian are given in Appendix 



IV. COMPARISON OF THE TWO APPROACHES 



We shall now compare the two models, the one using the effective Hamiltonian (fT3"]) derived from the master 
equation, and the phenomenological one using the Hamiltonian ([T5]) . We shall consider, as in [lfj, two pulses of the 
form 
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(16a) 



&2(t) = j^=sech (Jj; 1 



7T / t 

- tanh hi 

4 I T 



(16b) 



where T determines the widths of the pulses while a their intensities, so that aT gives a measure of each pulse area. 
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FIG. 2: (Color online) Asymptotic values of population P3 for the phenomenological and effective (perfectly superimposed) 
model for the intuitive sequence, as functions of T (in units of T^ 1 ), for aT = 10, AT = 1. 



A. Weak damping (IT < 1) 



Intuitive pulse sequence — For the intuitive sequence (i.e., when tt p = il± and fl s — ^2), the effective and 
phenomenological models deliver essentially the same results. This feature derives from the fact that the adiabatically 
transferred population is that of state |— ). In the weak-damping limit, the dominant source of decay is the population 
decay of |— ) itself due to the term (3,3) in both the Hamiltonians H c h and H p \ lcn , thereby diminishing the effects of 
the off-diagonal terms and in particular of ^Tsm2ip in (1,3) and (3,1), which are the only difference between the two 
models. Hence the decay of the hnal state, with a good approximation, is given by 



P.; 



exp 



2r 



sin 2 tp(t)dt 



(17) 



Figure[3]shows the asymptotic value of P3 as a function of T for both models, and shows the nearly perfect agreement 
of the two predictions, in particular the exponential dependence on T. 

Counterintuitive pulse sequence — For the counterintuitive sequence (f2 s = fii and fl p = Q.2), still under the 
weak damping conditions, we can study the evolution through the adiabatic elimination, assuming that the states 
|+) and |— ) are not very much populated during the process. We decompose the state vector in the adiabatic basis 
as = a+ |+) + ao |0) + a_ |— ), and the Schrodinger equation becomes a set of linear differential equations for the 
probability amplitudes, 



dt 



a + 




a+ 


a 


= H 


a 


a_ 







(18) 



with H = H c f[ or H = -ffphen, depending on the model considered. By setting d+ ~ d_ ~ 0, we find a + and a_ as 
functions of ao, and substitute these expressions in the equation for do, which assumes the form 



do(t) = -A{t)a {t), 

where A{t) has one of the following expressions depending on the model one is considering: 



(19) 



A»~(t) = ^ 2 (r + 2zOq cot 2ip) 

(iClo cot <p + r cos 2 ip) {—• iQo tan tp + T sin" tp) + ip 2 — |T 2 sin 2ip 



for the phenomenological model, and 



AeS = 2 [r(cos 4 tp + sin 4 y) + 2z» cot 2<p] 

(ifl cot ip + r cos 2 ip)(— iflo tan <p + V sin 2 1^3) + (p 1 



for the effective model. 
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FIG. 3: (Color online) Post-pulse values of population P3 for the phenomenological (solid red line) and effective (dashed blue 
line) model for the counterintuitive pulse sequence, as functions of T (in units of T -1 ) for weak damping only (a) and in a wide 
range (b). Relevant parameters: aT = 10, AT = 1. 



Equation (|TT)| has the very simple solution 

ao(i) = ao(— 00) exp 
which corresponds to the dark-state population 

P (t) = exp 



A{t')dt' 



For the phenomenological model, according to [10j, at the first order in the parameter T, we find 



3 (00) = exp 



-2r 



n 2 + ^ 



df 



while for the effective model the result is 

P 3 off M = exp 



;°° & (sinV + co3 4 y) dt ; 



\2 



(22) 



(23) 



(24) 



(25) 



where we have taken into account the fact that for the counterintuitive sequence, we have Po(~°°) = -Pi(~ 00 ) = 1 
and -Po(°°) = ^3(00). 

Since sin 4 cp + cos 4 cp < 1, it is evident that P| ff (oo) > P£ en (oo). Numerical results agree with this prediction, as 
shown in Fig. [3] 

It can be shown that as a increases the agreement between the results obtained in the adiabatic approximation 
and the numerical results improves. Indeed, generally speaking, higher a means higher Rabi frequencies between the 
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states, which put the system more under adiabatic conditions [7j. It is also worth nothing that the higher a the closer 
to unity are the values of P 3 for r = 0. Again, this depends on the validity of the adiabatic approximation, according 
to which the levels |+) and |— ) are not populated, so that the population transfer by |0) is more and more efficient, 
without losses of population to the other two levels. 



B. Strong damping: numerical simulations 

Figure [3}d shows the complete dependence of Pg ff (oo) and P|" hon (oo) on T and, in particular, the strong damping 
limit. While P3 >hen (cxD) decays as T increases, it is well visible that P| ff (oo) reaches higher values and in particular 
approaches 1 for high values of T. Such a behavior can be explained in terms of generalized quantum Zeno effect 
[2J, [25j|, meaning that the strong decay produces a separation of the Hilbert space into Zeno subspaces. If F is much 
larger than other quantities (f2o, ip, 9), we can split the Hamiltonian H c g as a sum of the unperturbed one, which 
contains only terms proportional to T, and a perturbation, that is 



H eS = 



—iT cos 2 ip 






fi[) COt if 

—id simp 

—iip 







— iT sin ip 

) sin ip iip 

—id cos ip 
'cos ip — fiotany? 



(26) 



The three eigenvalues of the unperturbed (first) part of the Hamiltonian are A = —iT cos 2 ip, —iT sir? ip and 0, cor- 
responding to the eigcnstates |+), |— ), |0). Since T is very high, for ip 7^ 7r/4 the three eigenstates correspond to 
very different eigenvalues and the presence of the perturbation only slightly changes the eigenvalues and eigenstates 
of the total Hamiltonian, so that the three states turn out to be essentially uncoupled throughout the process. On 
the contrary, when ip — the states |0) and |— ) could be coupled by the terms =pit ?cos ip. However, it is easy to see 
that when ip = 0, 9(t) is a constant function and 9 — 0. Moreover, when ip = ir/4 the states |+) and |— ) are coupled 
by the terms Hip, but the relevant transitions are unimportant since the only state with nonzero population is |0). 
Therefore, there are no transitions between the adiabatic states, and |0) can transfer the population from level |1) to 
level |3) without decay since the corresponding eigenvalue is essentially zero. 
In the phenomenological model, we have 



H, 



phen 



-iT cos 2 ip |rsin2<^ 



§rsin2<^ -irsin 2 <^ 

f^o cot ip iQ sin ip iip 
—i9 simp —iOcosip 

—iip iQ cos ip — flot&nip 



(27) 



and the eigenvalues of the unperturbed (first) part of the Hamiltonian arc A = 0, and — iT, corresponding to the 
eigenstates |0i) = |0), IO2) = simp |+) + cos ip | — ) and \—iT) = cos tp |+) — simp | — ), respectively. Since T is very large, 
the perturbation does not couple the state |— iT) to the doublet (the other two states in the degenerate subspace), 
but the states of the doublet can be coupled by the perturbation. In fact, the restriction of the perturbation to the 
doublet is 



doublet 



-it 
id 



(28) 



and since 9 in the counterintuitive sequence varies from to 7r/2, there is complete population inversion in the doublet . 
Therefore, since |0i) (—00) = |0) (—00) = |1), the final state of the system is | O2) (00) = |— ) (00) = |1), and there is 
no population transfer to level |3). Summarizing, in the strong coupling for a counterintuitive pulse sequence, the 
effective model predicts a complete population transfer whereas the phenomenological model predicts essentially no 
transfer. 

Let us now turn to the intuitive pulse sequence. We have already seen that for small values of T the population 
transfer from |1) to |3) is very low. This feature, which is common to the two models, relies on very different physical 
mechanisms. Let us start considering that in the intuitive sequence, for the phenomenological model, at t = —00 we 



8 



(a) 






J » 







6 







4 







2 





-4-2 2 4 



t 

(b) 

1 r 

0.8 
0.6 - 
0.4 - 
0.2 - 

^ ^ ^ ^ 

-4-2 2 4 

t 

FIG. 4: (Color online) Evolution of the populations in the strong damping limit, according to the effective model in the intuitive 
sequence (a) and in the phenomenological model in the counterintuitive sequence (b): Pi as a solid (red) line, P2 as a dotted 
(green) line and P3 as a dashed (blue) line. Populations P2 and P3 are essentially zero everywhere. Time is measured in units 
of T. Relevant parameters are oT = 10, FT = 500, AT = 1. 

have |1) = |— ) = [O2) , which, according to the previous analysis, undergoes a complete transition toward state |0i), 
so that the final state is |0i) (00) = |0) (00) = and there is no transfer to state |3). Figurc@Ji shows this behavior. 
In the effective model, the system starts from |— ) (—00) = |1) and there are transitions between states |+) and |— ) 
around ip — 7r/4; the entire subspacc spanned by these two states decays. Therefore, a complete loss of probability of 
the system characterizes the dynamics in this case. Figure [4b illustrates this feature. 

V. CONCLUSIONS 

In this paper we have shown that the phenomenological and effective models for a laser-driven three-state A-system 
under dissipative dynamics give different descriptions. 
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We have shown that the two relevant non-Hcrmitian Hamiltonian models differ, in the adiabatic basis, in off- 
diagonal terms which couple levels 1 and 3 and which are proportional to the decay constant T. This suggests that 
the discrepancy between predictions coming from the two models may increase as the decay rate increases. Both 
analytical and numerical results confirm this insight. Indeed, for weak damping the predictions from both models are 
essentially the same, while in the strong damping limit the difference is evident. Specifically, in the strong damping 
limit, for the counterintuitive pulse sequence, there is a significant difference in the values of the post-pulse population 
of the target state: according to the phenomenological model there should be no transfer from state 1 to state 3, 
while the effective model predicts almost complete population transfer as T increases. We note that, even though for 
the intuitive pulse sequence the predictions for the post-pulse population of level 3 are the same for both models, the 
physical mechanisms that lead to this result are very different for large T. In fact, in the phenomenological model the 
population is kept in state 1, while in the effective model there is no final population in state 3 because all the states 
have undergone a total decay. 

The most important result in this paper is the prediction of a complete population transfer for the counterintuitive 
sequence in the strong damping limit; this result is of potential interest in applications of STIRAP schemes in the 
manipulation of quantum states. 
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Appendix A: Derivation of the Effective Model 

In this appendix we derive the effective Hamiltonian model in (|13|) from the master equation (J7|). 
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Substituting ((9|) into (J7J, and assuming zero temperature (7's = 0) one obtains the following set of rate equations 
in the adiabatic basis: 



Poo = 
P++ 
P— 
P+o = 

P+- 



1 sin p p+o — v sin p p 0+ — 6 cos tp p_o — " c °s p po- , 
7+ P++ + <P P-+ + (pp+-+0smip p + + sin p p +Q , 
7-/? — — tp p-+ — tp p+~ + cos p p 0+ + 9 cos tp p +0 



7+ 



p +0 + 8 sin tp poo + pp-o — sin p p ++ — 9 cos p p + . 



Po- 



-IWo- 



7+ +7- 



Po- 



p+- + 6 sin p po- + pp — + 9 cos p p +0 - tp, p++ , 

1 sin p> p + - - 9 cos 4> p — — tp p Q+ + 9 cos p p o — P P++ , 



(Ala) 
(Alb) 
(Ale) 
(Aid) 

(Ale) 

(Alt) 



and the Hermitian conjugates of the last three equations. Equations involving level 4 are not shown, since they don't 
play any role in the effective model. As a consequence the effective model does not conserve the total probability. It 
is straightforward to see that substituting the Hamiltonian model in (TT5)) into (|T^1) , one obtains exactly the same rate 
equations. Therefore we conclude that H e g represents the effective Hamiltonian of the system. 



Appendix B: Strong damping in the bare basis 

In this appendix we show how to treat the strong damping dynamics in the bare basis {|1) , |2) , |3)}. In ref. [To| 
the starting point of the treatment is to give the phenomenological Hamiltonian: 



Tjbare 
phen 



n p (t) 

fip(t) A - iT n s (t) 







fi,(t) 







(Bl) 



which, transformed to the adiabatic basis {|+) , |0) , |— )}, gives Eq. ((T5|) . Equation (|B1[) clearly indicates that in the 
strong damping limit state \2) is well separated from the other states, so that the coupling scheme does not allow to 
transfer population from state |1) to state |3). 
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In the master equation approach, the effective Hamiltonian in the bare basis assumes a much more complicated 
form which does not allow to separate the three bare states. Indeed, since the effective Hamiltonian is related to the 
phenomcnological one by the relation 



H, 



eff 



H, 



phen 



2 Y 



1 



1 



(B2) 



we find that the inverse transformation from the adiabatic to the bare basis gives: 



r 

P hcn"*2 sin 2^ 



2 sin tp cos ip sin cos 2p sin 6 2 sin <p cos ip sin 8 cos ( 
cos2ips'm9 — sm2p cos2</?cos# 

2 sin ip cos ip sin cos 6 cos 2ip cos 9 2 sin ip cos p cos 2 6 



(B3) 



which clearly shows the impossibility of separating state |2) from the other ones. Therefore the only basis in which 
the strong-damping dynamics can be easily explained is the adiabatic one. 
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